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Abstract 

This paper introduces an extension of the Markov switching ARCH model 
where the volatility in each state is a convex combination of two different ARCH 
components with time varying weights with different volatilities. The asymp- 
totic behavior of the second moment is investigated and an appropriate up- 
per bound for it is evaluated. The estimation of the parameters by using the 
Bayesian method via Gibbs sampling algorithm is studied. We propose a dy- 
namic programming algorithm for the forecasting. Finally we illustrate the 
efficiency of the model by simulation and forecasting the volatility. We show 
that this model provide much better forecast of the volatility than the Markov 
switching ARCH model. 
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1 Introduction 

In the past three decades, there has been a growing interest in using non linear time 
series models in finance and economy. For financial time series, the ARCH model 
and GARCH model , introduced by Engle pT] and Bollerslev [7], are surely the 
most popular class of volatility models. Although these models have been applied 
extensively in the modeling of financial time series, but the dynamic streucture 
of volatility can not be captured passably by such models. For more consistent 
volatility modelling, the models by time varying parameters are introduced. One 
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class of such models is that of smooth transition GARCH models that presented 
by Gonzalez-Rivera [15], (see also Hagerud [20] and Medeiros and Veiga [26]). An- 
other class is that of Markov switching models. These models are obtained by by 
Merging (G)ARCH model with a Markov process, where each state of the Markov 
model allows a different (G)ARCH behavior. These models are introduced by Cai 
[8] and Hamilton and Susmel [19]. This feature extends the dynamic formulation 
of the model and potentially enables improving forecasts of the volatility pQ. Gray 
[16j . Klaassen [22], Haas, Mittnik and Paolella [18] proposed different variants of 
Markov-Switching GARCH models. See also further studies, Abramson and Cohen 
[JJ, Alexander and Lazar [2] and Bauwens et al. [6]. The component GARCH mod- 
els, introduced first by Ding and Granger [10], are also a generalization of constant 
parameter GARCH model. These models have been widely applied in modeling the 
financial time series (e.g. [12], |25] and [J3]). In the structure of component GARCH 
model ([10]), two different ARCH component contribute to the overall conditional 
variance at time t. One component has the high volatility (integrated variance com- 
ponent) and the other component has the low volatility. A generalization of the 
component GARCH model of Ding and Granger is the weighted GARCH model 
that is peoposed by Bauwens and Storti [5]. In this model the weights of GARCH 
components are the function of lagged values of the conditional standard deviation 
or squared past observations. 

In this paper we consider a Markov switching model that the volatility of each 
state is a convex combination of two ARCH regimes with time varying coefficients 
which is in effect of previous observation. This model has the potential to model 
the effect of more complicated resources which are in effect of some volatility com- 
ponents and the share of these components could change in time. We consider 
different weight functions for each state that allow volatility in each state to react 
differently to the shocks of equal size. As using all past observations for forecasting 
could increase the complexity of the model, we reduce the volume of calculations by 
proposing a dynamic programming algorithm. We derive necessary and sufficient 
conditions for stability and obtain an upper bound for the limit of the second mo- 
ment by using the method of Abramson and Cohen [1] and Medeiros [26]. For the 
estimation of the parameters we use the Bayesian inference via the Gibbs sampling. 
We compare the performance of our model to Markov switching ARCH model. The 
Markov switching component ARCH model can forecast the conditional variance 
much better than MS- ARCH model. 

The paper is organized as follows: in section 2 we introduce the smooth transition 
Markov switching ARCH model. Section 3 investigates the statistical properties of 
the model. Section 4 is devoted to estimation of the parameters of the model. 
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Section 5 dedicated to the analyzing of the efficiency of the proposed model through 
simulation and the comparison of the forecast errors with the MS-ARCH model. 
Section 6 concludes. 



2 Markov switching Component ARCH model 

The Markov switching component ARCH model, MS-CARCH, for time series {yt} 
is defined as 



Vt = £t\/ H t,z t > H t,z t = w t ,z t hi, t ,z t + (1 - w t ,z t )h2,t,Zt, (2.1) 

where {et} are iid standard normal variables, {Zt} is an irreducible and aperiodic 
Markov chain on finite state space E = {1,2, ••• ,K} with transition probability 
matrix P = \\pij\\KxK, where -pij = v{ z t = j\ z t-i = i), hj G {!,••• ,K}, and 
stationary probability measure ir = (m, • • • ,ttk)'- Also 

hi, t ,z t = a 0tZt + ai,z t yl-x, h 2 ,t,z t , = b ,z t + h,ZtVt-i> ( 2 - 2 ) 

and each of the weights (wti,i = 1, • • • , K) is a function of past observation as 

1 -exp(-7j|j/ f _i|) 

wt,i = ~, : - f i rr 7i > °> ( 2 - 3 ) 

l + exp(-7i|y t _i|) 

which is bounded , < u^t-i < 1. The parameter 7, is called the slope parameters, 
that explains the speed of transition from one component to the other one: the 
higher 7^, the faster the transition. Since 7, > 0, when the absolute value of yt-i 
increases, the impact of h± i t is increases and consequently the effect of /i2,i,t decreses 
and vice versa. For this reason we consider the first ARCH component in each state 
with the high volatility and the second component with low volatility. So when 7 
tending to zero or infinity and the MS-CARCH model tends to MS- ARCH model. 

It is assumed that {et} and {Zt} are independent. Sufficient conditions to guar- 
antee strictly positive conditional variance are ao,j, 6o,i to be positive and a\^, a>2 t i, 61^, 62,. 
being nonnegative. 

Let It be the observation set up to time t. The conditional density function of 
yt given past information is obtained as follows: 
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A 



f(y t \l t -i) = ^2 f(Vt> Zt = ifa-i) 
i=i 

K 

s ^ j p{Z t = i\Xt-i)f(y t \Zt-i,Zt 



i=l 
K 



E-fV(^=) (2.4) 



in which otf 1 = p(Zt = i\Xt—i) (that is obtained in next section), and <fi(.) is the 
probability density function of the standard normal distribution. 



3 Statistical Properties of the model 

In this section, the statistical properties of the MS-CARCH model are investigated 
and the conditional variance of the process is obtained. We show that the model, 
under some conditions on coefficients and transition probabilities , is asymptotically 
stable in the second moment. An appropriate upper bound for the limiting value of 
the second moment is obtained. 



3.1 Forecasting 

The forecasting volatility (conditional variance) of MS-CARCH model is given by 

K K 

Var{Y t \Jt-i) =Y, a T H t,i = X^fK^W + (! - ™t,i)h 2 ,t,i) (3.5) 
i=i i=i 

as Hk t t is the conditional variance of k-th state. This relation shows that the con- 
ditional variance of this model is affected by the changes in states, the volatility of 
components and the weight functions in each state. 

At each time t, ctf* (in equation (|2.4p . (|3.5|) ) can be obtained from a dynamic pro- 
gramming method based on forward recursion algorithm, proposed in remark (|3,ip . 



Remark 3.1 The value of is obtained recursively by 



(t) = Em=i f(yt-i\Z t -i = m,Zt-2)p{Zt-\ = m\l t _ 2 )p 
E ro =i f(yt-l\Zt-i = m,Xt-2)p(Z t -i = m\T t -2 
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Proof 3.1 As the hidden variables {Zt}t>i have Markov structure in MS-CARCH 
model, so 

K 

M 



where 



ay =p(Z t = j\l t -i) = ^ P(Z t = j, Z t -i = m\l t -i) 

m=l 

K 

= ^2 P( z t = j\ z t-i = w,It_i)p(Z t _i = m\It-x) 



m=l 
K 



p{Z t = j\Z t -i = m)p(Z t -\ = rn\It-i] 

m=l 

Em=l fi X t-i, Z t -i = m)p mJ 



E m =i/(2t-i»^t-l = m) 



Em=i f(yt-i\Z t -i = m,l t -2)p(Z t -i = m\Zt- 2 )p 



Em=i /(yt-ll^t-i = rn,l t -2)p(Zt-l = m\l t - 2 ) 



f{y t -\\Z t -i = ?n,X t _ 2 ) = 



t—l,m 



(3.7) 



3.2 Stability 

In this subsection, we investigate the stability of second moment of MS-CARCH 
model. Indeed we are looking for an upper bound for the second moment of our 
model. The second moment of the model can be calculated as: 

E(y 2 t ) = E(H tiZt ) = E Zt [E t ^(H tjZt \z t )] 



K 

= Y,Kz t Et-i(H t , Zt \z t ). (3.8) 

zt=l 

Et(-) denotes the expectation with respect to the information up to time t. Also 
for summarization, we shall use E(-\zt) and p(-\zt) to represent E(-\Zt = zt) and 
P(-\Zt = zt), respectively, where zt is the realization of the state at time t. We 
investigate E t -i[Ht : Zt\ z t] as follows: 



5 



Et^Ht^zt) =E t -i{[w t , zt {aQ jZt + ai,^ t y t 2 _i) + (1 - w t , Zt )(b , Zt + &i,*|/t_i)]|*t) 
= b , Zt +h,zt E t-i[yt-i\ z t\ + ( a o,2* - b 0>Zt )Et-i[wt,z t \zt] 

^-V-^ » v ' ' v ' 

/ 77" 777 



IV 

(3.9) 



The relation (II) in (|3.9p can be interpreted as follows: 

A 



Z t -i=l JS Zt-l 



t-i\zt, Zt-i)p(zt-i\zt)dZt- 



K 



= ^2 p( z t-i\ z t) E t-i[yt-i\ z t-i,zt], (3.io) 

2t-l=l 

where S , 2 t _ 1 is the support of It-i = (yi, • • ■ , Vt-i)- Since the expected value of 
conditional on the present state is independent of any future state, so 

E t -i[yt-i\zt-i,zt] = £tf-i[j/t_i|zt-i]. (3.11) 

Also using the tower property of the conditional expectation, that is E[E(Y\X, Z)\X] 
E(Y\X) [see Grimmett and Stirzaker (2001, p. 69)], we have 

E t -i[yt-i\z t ~i} = E t _ 2 [Et-i(yt-i\Zt-2,zt-i)\zt-i} 

= Et^Ht^z^zt-i]- (3-12) 

The calculation of Et-i[wt, Zt \ z t\ and Et-i[wt :Zt y^_i\zt] is a problem that can not 
be easily done, For this reason we will try to find an upper bound for them. 

Upper bound to III. As < w^i < 1, so an upper bound for the relation III in 
(13. 9|) is obtained by 



(ao,z t ~ bo,z t )Et-i[wt tZt \zt\ < {oq, h - b ,z t ) < oo. (3.13) 
Upper bound to IV. Let < M < oo be a constant, so 

Et-iiwt^Vt-ilzt] =E t - 1 [wt, Zt yt-iI\y t _ 1 \ < Ai\zt} 

+ E t -.i[w t j lt yi-iI\ vt _ 1 \>M\zt] 
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in which 

_ j 1 if x < a 
x<a \ otherwise. 

As by (f!T3|) . < w t)Zt < 1 and so 
also 

Et-l[wt,z t 1/t.l I \yt-i\>M\ z t] = / I/iiKzJpft-il 2 *)^-! 

J 'Sx t _ 2 ,y t -i<-M 

+ / y^iK^M^-iN^-i, 

JSz t _ 2 ,yt-i>M 

by (JZ3J, 

lim io t , = 0, lim w ttZt = 0, (3-14) 

!/t-l-H-oo j/t— i— oo 

therefore according to the definition of limit at infinity, for a small number 5 > 0, 
there will exist a finite constant M > such that if yt-i > jW, l^zj < 5 and if 
S/t— l < -M, [io t % | < 5. Hence 



^-l^tJ/t-i^ii/t-il^k] < <m yt-iP(Zt-i\zt)dZt-i 

+ <5 / yl_ x p(It-x\zt)dZt-\- 

J St. „,j»-i>M 



Since the distribution of the {et} is symmetric, then 
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/ y 2 _ x p(I t -i\zt)dXt-i <5 / y\_ x p{&-\\zt)$Lt-\ 

J Si t _ 2 ,y t -i<-M JSx t _ 2 ,-oo<y i _i<0 



and 

5 / y$_ 1 p(Z t -. x \z t )dZt-.i <5 / yf^pilt-ilz^dlt-i 

J S Xt _ 2 ,yt-i>M J S It _ 2 ,0<yt-i<oo 



6- 
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Therefor 

^-iK*!/t-iN <M 2 + SEt-i^M. 
By replacing the obtained upper bounds and relations (I3.10p - (l3.12p in (13, 9D . the 
upper bound for E t -i{H t .z t \zt) is aquired as: 

E t -x{H tj z t \zt) < ao, Zt + (ai, 2t - & Mt )M 2 

K 

+ ^ [&i,*t + (ai,et - &i^)^pfe-i|^)#t-2[#t,z t _i|^-i], (3-15) 

zt-l=l 

in which by Bayes' rule 

p(z t -i\z t ) = 1 ^ t± {Pz t - i z t }, 
where P is the transition probability matrix. Let 

n = [oo,i + (ai,i - 6i,i)M 2 , • • • , a ,fi: + Kx - &i,*r)M 2 )]', (3.16) 
be a vector with K component, C be a K-by-K matrix with elements 

{C jk } = [hj + (aij - &ij)<J]— {P fei }, (3.17) 

and 

A t = [£^i(ff ti i|Z t = l) J .-- ,£ t _iOfV|^ = (3.18) 
be a K-by-1 vector. 

Hence by (I3.16p - (13.18|) we have the following recursive inequality, 

A t <n + CA t _i, *>0. (3.19) 
with some initial conditions A_i. 

Suppose p(A) denotes the spectral radius of a matrix A, then we have the following 
theorem. 

Theorem 3.1 Let {Y t }^ follows the MS-CARCH model, defined by IE7]) - IO|) . 

the process is asymptotically stable in variance and lim^oo E(Y t 2 ) < ir'(I — C)fi if 
and only if p(C) < 1. 

Proof 3.2 By recursive inequality i3.19\) , 
t-1 

A t <tlJ2 ci + ctA o : = B t 

i=0 

(3.20) 
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By the matrix convergence theorem 123)/ , a necessary and sufficient condition for 
the convergence of Bt when t — > oo is p(C) < 1 ( the value of 5 can be considered 
small enough to be negligible). Under this condition, C* converges to zero as t goes 
to oo and Yli=o converges to (I — C)" 1 provided that matrix (I — C) is invertible. 
Soifp(C)<l, 

lim A t < (I - C) _1 « 

t— >oo 

and by \3. 8\) we attain the upper bound for the asymptotic behaviour of unconditional 
variance, 

Hm t ^ooE(yf) < tt'(I - C)n. 



4 Estimation 

In this section we describe the estimation of the parameters of the MS-CARCH 
model. We consider Bayesian MCMC method using Gibbs algorithm by following 
methods of sampling of a hidden Markov process ([9] and |21j). MS-GARCH model 
and weighted GARCH model ([5] and [6] ) for estimation of parameters. 

Let Y t = - , yt) and Z t = {z\,--- ,zt). For the case of two states, the 

transition probabilities are n = (rju, 7712, %i) V22) and the parameters of the model 
are = (0i,0 2 ), where 6 k = (a ok ,bok,ai k ,bi k ,^ k ) for k = 1,2. 

The purpose of Bayesian inference is to simulate from the distributions of the 
parameters and the state variables given the observations. As Z = (z\, ■ ■ ■ , zt) and 
Y = (yi, ■ ■ ■ , yr) the posterior density of our model is: 

p(0, 17, Z\Y) oc p(9, f])p(Z\e, r,)f(Y\9, r,, Z), (4.21) 

in which p(6,rj) is the prior of the parameters. The conditional probability mass 
function of Z given the (9, rj) is independent of 6, so 

p{Z\9,n) =p(Z\r)v< h mi) 

T 

= np{zt+i\zt, Voo,Vii) 
t=\ 

= Poo°°(l -Poo) n >? 1 11 (l -Pii) ni °, (4.22) 
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where = #{zt = j\zt-\ = i}. The conditional density function of Y given the 
realization of Z and the parameters is factorized in the following way: 



T 



f(Y\ V ,d,Z) = Hf(y t \d,z t = k,Y t _ l ), k = 1,2, (4.23) 



t=l 

where the one step ahead predictive densities are: 

f(y t \d,z t = k.Yt-i) = -=L=exp(--^-). (4.24) 

^]l'Ktit,k M t,k 

Since the posterior density (|4.2ip is not standard we can not sample it in a 
straightforward manner. Gibbs sampling of Gelfand and Smith [T3] is a repetitive 
algorithm to sample consecutively from the posterior distribution. Under regularity 
conditions, the simulated distribution converges to the posterior distribution, (see 
e.g Robert and Casella [27]). The blocks of parameters are 6, i] and the realizations 
of Z. 

A brief description of the Gibbs algorithm: Let use the superscript (r) on Z, 9 and 
r\ to denote the estimators of Zt), rj, and 9 at the r-th iteration of the algorithm. 
Each iteration of the algorithm consist of three steps: 

(i) Drawing an estimator random sample of the state variable Z^ given , rj^ r ~ 1 \ #( r ~ 1 ). 

(ii) Drawing a random sample of the transition probabilities r/ r ) given 
(hi) Drawing a random sample of the 0^ given Z r and r/ r ). 

These steps are repeated until the convergency is obtained. In what follows 
sampling of each block are explained. 



4.1 Sampling z t 

The purpose of this step is to obtain the sample of p(zt\r],8,Y t ) that is performed 
by Chib[9], (see also [H]). Suppose p(zi\r), 9, Yq, ) be the stationary distribution of 
the chain, 

P {z t \r 1 ,9,Y t )(xf(y t \9,z t = k,Y t _ l )p{z t \r 1 ,9,Y t _ l ), (4.25) 

where the predictive density f(y t \9,zt = /c, Yj_i) is calculated by the relation (|4.24p 
and by the law of total probability p(zt\r],9,Y t -i) is given by: 

K 

p(z t \ V ,9,Y t _ 1 )= K^-iM^t-iK-m- ( 4 - 2 6) 

Zt-l=l 
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Given the filter probabilities (p(zt\r],9,Y t )), we run a backward algorithm, start- 
ing from t = T that zt is derived from p(zt\ij, 0, Y). For t = T —1, ■ ■ ■ , the sample 
is derived from p(z t \z t+ i, ■ ■ ■ , zt, 0, 77, Y),which is obtained by 

p{z t \z t+1 ,- ■ ■ ,z Tl 0,r],Y) oc p(z t \r],e,Yt)r] ZuZt+1 . 

To derive Zt from p(zt\-) = p Zt is by sampling from the conditional probabilities 
qj = p(Z t = j\Z t > j, .) which are given by 

p(Z t =j\Z t >j,.) = -^—. 

After generating a uniform (0,1) number U, if U < qj then zt = J, otherwise increase 
j to j + 1 and generate another uniform (0,1) and compare it by c^+i- 

4.2 Sampling rj 

This stage is devoted to sample rj = (r/u , 7722) from the posterior probability p(rj\9, Y t , Z t ) 
that is independent of Yj, 0. We consider independent beta prior density for each of 
rin anci V22- F° r example, 

pfaii|Z t ) oc p( m i)p(Z t \rm) = V ^ +nil - 1 (l - r) n ) c ™ +n ™-\ 

where en and C12 are the parameters of Beta prior, is the number of transition 
from zt-i = i to zt = j. In the same way the sample of 7/22 is obtained. 

4.3 Sampling 9 

The posterior density of 9 given the prior p(9) is given by: 

T T 1 2 

p(0|y, Z, r,) oc p(0) /G/tl*. *t = *, ^-l) = P(0) U nr^- exp(--p-), (4.27) 

which is independent of 77. Since the conditional distribution of 9 does not have a 
closed-form (because for example p(aok\Y t , Z t , 9- aok ), in which 9_ aQk is the param- 
eter vector without aofc, contains Hk t i, which is also a function of aofc- Therefor it 
isn't a normal density.) using the Gibbs sampling in this situation may be compli- 
cated. The Griddy Gibbs algorithm, that introduced by Ritter and Tanner (1992), 
can be a solution of this problem. This method is very applicable in researches (for 
example [4] , [5] and [6]). 

Given samples at iteration r the Griddy Gibbs at iteration r + 1 proceeds as 
follows: 
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Figure 1: Simulated time series of MS-CARCH model. 
Table 1: Descriptive statistics for the simulated data (sample sizc=200) 



Mean 


Standard deviation 


Skewness 


Maximum 


Minimum 


Kurtosis 


.051 


.632 


.445 


2.893 


-1.833 


5.777 



1. Select a grid of points, such as ajj, a^, • • • , a^. Using (|4.27p . evaluate the condi- 
tional posterior density function k(aoi\z t ,Y t ,0- ao ) over the grid points to obtain the 
vector G k = {fa,-- - ,k G ). 

2. By a deterministic integration rule using the G points, compute G$ = (0, $2, • " " > 
with 

^ = &(a 01 ^ , z\ r) , Y t )da 0i , i = 2, ■ ■ ■ , G. (4.28) 

Ja 0i 

3. Simulate u ~ ?7(0, 3>g) and invert <l?(aoi|#^a 0i , -^ r \ Y" 4 ) by numerical interpola- 
tion to obtain a sample aoi +1 ' > f rom o-0i\G-a 0i ^ %t ■> Y t . 

4. Repeat steps 1-3 for other parameters. 

For the prior densities of all elements of 9, it can be can considered independent 
uniform densities over the finite intervals. 



5 Simulation results 

In this section we provide some simulation results of MS-CARCH model defined 
by equations f)2. l|) - ()2.3j) for two states. We simulate 200 sample from the following 
MS-CARCH model: 

y t = e t ^mj t , (5.29) 
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Figure 2: Squared observations of the simulated time series (blue), forecasts by MS- ARCH model (red) 
and forecasts by MS-CARCH model (green) 
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Figure 3: Forecast Errors of square of the observations in the MS- ARCH model (red) and in MS-CARCH 
model (green). 
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Table 2: Results of the Bayesian Estimation of the simulated MS- ARCH model 





T^TllP Trainee; 


1\/Tps} n 


Std Hpv 




0.10 


0.097 


0.049 


an 


0.40 


0.397 


0.164 


&01 


0.20 


0.199 


0.053 


hi 


0.90 


0.845 


0.190 


a-02 


0.10 


0.158 


0.081 


a\2 


0.05 


0.124 


0.067 


bo2 


0.30 


0.298 


0.110 


h2 


0.20 


0.200 


0.108 


7i 


5 


5.994 


.272 


72 


0.10 


0.106 


0.057 


Vn 


0.90 


0.892 


0.078 


r}22 


0.85 


0.709 


0.159 



where {et} is an iid sequence of standard normal variables, {Zt} is a Markov chain 
on finite state space E = {1,2} with transition probability matrix 



P 

and 



.90 .10 
.05 .85 



H tA =(1 - Wt>1 )(.l + AyU) + m,i(-2 + .9yU), w t>1 = ] GXp{ 

1 + exp(-5|y t -i|) 

1 - exp(-.l[y t _i| 



H t)2 =(1 - w tt2 )(.l + .Qby^) + w t>2 (.3 + .2^), w tj2 



1 + exp(-.l|y t _i|) 

(5.30) 



The first state implies a higher conditional variance than the second one and 
in each state, the first component has the lower volatility than another component. 
By theorem (1), the assumption for the exictence an upper bound for the second 
moment is checked. 

In table 1, we report summery statistics for simulated data and figure 1 shows 
the plot of the simulated time series. 



Using the Bayesian inference, we estimate the parameters of the MS-CARCH 
model. The prior density of each parameter is assumed to be uniform restricted 
over a finite interval (except for ijn and r]2 2 , since they are drawed from the beta 
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distribution). The number of iterations of the Gibbs sampler was 20000 and the 
initial 10000 draws were discarded. Table 2 demonstrates the performance of our 
estimation methods for the model. The results of these tables show that the stan- 
dard deviation are small enough in most cases. 

For clarifying the performane of MS-CARCH model toward MS-ARCH model, 
We compare the forecasting volatility (E(Y^\Tt~i)) of each model to the squared 
observations. The forecast error (the difference between the forecasting volatility 
and the squared observations) of our model is more smaller than the MS-ARCH 
model especially in picked points. 

6 Conclusion 

In this paper a generalization of the MS-ARCH model has been presented where the 
conditional variance in each state is a convex combination of two different ARCH 
components with time varying coefficients, one of the component with higher volatil- 
ity than other component. Our model can providemore better forecast of volatility 
toward MS-ARCH model. For the estimation of parameters we have applied the 
Bayesian estimation algorithm. We provide simple necessary and sufficient condi- 
tion for the existence of an upper bound for second moment. 

This work has the potential to be applied in the context of financial time se- 
ries. The empirical distribution of daily returns doesn't generally have a Gaussian 
distribution. They have fat tails densities (they are called leptokurtic) . One of 
the extending of this work is considering the fat tail densities instead of Gaussian 
distribution, that can cause better modeling of the financial time series. Also we 
can generalize this model by using of the GARCH structure instead of the ARCH 
structure for the better results. 
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